Aeroelastic instability of a flexible plate in a uniform flow 
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The flutter instability of a flexible plate immersed in an axial flow is experimentally investigated. 
This instability is similar to flag flutter and results from the competition between destabilizing 
pressure forces and stabilizing bending stiffness. The instability threshold is measured as a function 
of the experimental parameters (plate size and aspect ratio) and compared to different theoretical 
predictions. 
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The flutter of a flexible plate immersed in an axial flow 
is a canonical example of flow-induced vibrations. This 
instability can be experienced in everyday life when one 
observes a flag flapping in the wind. Because this phe- 
nomenon appears in many applications (paper industry, 
airfoil flutter, snoring), it has motivated a large litera- 
ture which has been recently reviewed [ij. This instabil- 
ity can be regarded as a competition between fluid forces 
and elasticity. Indeed, as soon as the plate experiences a 
small lateral deflection, a destabilizing pressure jump ap- 
pears across the plate, while the bending stiffness tends 
to bring the plate back to the planar state. 

This system can be studied by restricting the analy- 
sis to one-dimensional flutter modes as observed in most 
experiments. In this case, the plate motion obeys the 
Euler-Bernoulli beam equation with additional pressure 
forces which are calculated by assuming a potential flow. 
As a flrst approach, theplate can be supposed to be infi- 
nite in both directions [2| to calculate the fiow. In other 
theoretical studies, the plate length L (or chord) takes 
a finite value while two asymptotic limits are considered 
for its span H . If 7? ^ i, the fluid forces can be calcu- 
lated using the Lighthill's slender body theory [3J, Ij and 
\i H ^ L the problem can be treated as two-dimensional 
[a, la, lZI, y. In this latter case, the flow is entirely de- 
scribed by point-vortices which are distributed within the 
plate and possibly in its wake. It is known from airfoil 
theory that this problem does not admit a unique solu- 
tion (intrinsically because the Laplace equation has to 
be solved on an open domain). This non-unicity can be 
treated by two approaches [5(]: one can either consider a 
zero total circulation around the plate and thus no vortic- 
ity in the wake, or the circulation can be prescribed by 
the Kutta condition applied at the trailing edge of the 
plate [SI] . This latter condition imposes advected vortices 
in the wake which makes this model more physical for 
this flag-type instability. Note that this shed vorticity 
can equivalently be regarded as the consequence of the 
Kelvin's circulation theorem. 

Using these two-dimensional flow models, a stability 
analysis can be carried out that allows to predict the 
critical velocity for plate flutter. These theories have 
been compiled by Watanabe et al. [6| and compared with 



existing experimental data [3|, |7| . Remarkably all these 
theories predict approximately the same critical velocities 
but strongly underestimate the measured thresholds. In 
other words, the plate appears systematically more sta- 
ble than predicted by a two-dimensional approximation. 
This discrepancy has motivated the recent study of Eloy 
et al. jlQj in which the finite plate span is explicitly taken 
into account in the analysis. The present paper aims at 
comparing their predictions with experimental measure- 
ments in which the effect of the plate aspect ratio H/L 
is extensively investigated. 

As shown in Fig.[l] we consider a flexible plate of span 
H and length L, lying in the vertical plane (Oxy) and im- 
mersed in an axial flow of velocity U. Its flexural rigidity 
is given by D = Eh'^ /12{1 — i/^), where E is the Young's 
modulus, h the plate thickness and i/ its Poisson's ratio. 
For small lateral deflections z{x,t), the plate motion is 
driven by the linearized Euler-Bernoulli equation 



mdfz + Dd^z + {[p]) = 0, 



(1) 



where m is the mass per unit surface of the plate, 
[p]{x,y,t) is the pressure jump across the plate and the 
notation (.) stands for the average along the span H. 
Equation ([1]) is valid in the limit of an inviscid fluid 
and for a material of negligible visco-elastic damping. 
In addition the deflection z satisfies clamped boundary 
conditions at the leading edge: z{0) = dxz{0) — and 
free boundary conditions at the trailing edge: d'^z{L) = 
dlz{L) = Q. 

Assuming an inviscid flow, the perturbation velocity 
can be fully described by the vorticity distribution in the 
flow. Vortex-lines are located in the plate and in its wake 
as sketched on Fig. [TJb). The pressure jump across the 
plate only depends on the y-component of the vorticity 
j{x,y,t) through the unsteady Bernoulli equation 



d^[p]^pu{dt + ud,)j, 



(2) 



where p is the fluid density. As seen from this equation 
a vorticity distribution of the form jo{y) exp[i(fca; — tot)] 
with wave velocity co/k — U is compatible with a zero 
pressure jump and is indeed the distribution of vorticity 
in the wake. 
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FIG. 1: Sketch of the experimental setup (a); and plate sub- 
ject to the one-dimensional deflection z{x,t) (b). 



In order to couple the vorticity distribution 7 to the 
plate deflection z, a kinematic condition is enforced on 
the plate surface which can be written as 



{^mK{x - Ode = {dt + Ud,)z = w{x, t), (3) 



where ^ is a dummy variable in place of x, w is the normal 
plate velocity in the fluid reference frame and the letter 
C across the integral sign indicates that the Cauchy prin- 
cipal value should be taken. Since this inverse problem 
can have non-unique solutions, the Kutta condition is 
also used at the traihng edge i.e. \p] (L) = 0. The kernel 
K in Eq. ^ expresses the influence of a vertical vorticity 
line located in ^ on the z-component of the fluid veloc- 
ity in X. For an infinite span, the azimuthal velocity of 
a point-vortex in 2D yields the kernel K{X) = 1/2ttX 
used in 2D models [5|. For an asymptotically small H 
the Lighthill's slender body theory is equivalent to tak- 
ing K{X) = sgn{X)/TiH. In Eloy et al. 0, we have 
shown that an approximate kernel can be used for any H 
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sgn(X), (4) 



which is constructed to have the correct asymptotic be- 
havior for small and large H. 

To address the stability of this fluid-structure cou- 
pling, the full linearity of this problem is used and 
Galerkin decomposition can be assumed for the deflec- 
tion: z(x^t) = ^„ 2;„(a;) exp(itji). For a given Galerkin 
mode, one has then to invert Eq. ([3]) and inject the solu- 
tion in Eq. ^ to find the corresponding pressure jump 
([p„]) exp(iwi). The partial differential equation ([T]) is 
then reduced to an eigenvalue problem for the complex 
frequency lo whose eigenmodes correspond to the insta- 
bility modes. For small flow velocity, all modes are stable 
(i.e. their complex frequencies w have a positive imagi- 
nary part). When U is larger than a critical flow velocity 
C/c, one of the instability modes eventually becomes un- 
stable. 



Using L and LjXJ as characteristic length and time, 
the system parameters are reduced to three dimensionless 
numbers: the reduced velocity C/*, the mass ratio M* and 
the aspect ratio H* given by 



U* = LU^Jm/D, M* = pL/m, H* = H/L. (5) 

In this dimensionless form, the present theory (see 
Ref. [lOl for details) allows to predict the critical veloc- 
ity U*{M* ,H*) and the dimensionless mode frequencies 
uj* — lliL/U. Note however that the pressure singular- 
ity at the leading edge has been artificially suppressed 
in this analysis in such a way that that, in the limit of 
infinite H* , the results of the 2D theory of Kornecki et 
al. [5i] are not exactly recovered (their results have been 
recalculated in the present paper with better computer 
accuracy to allow comparisons). 

As sketched in Fig. [Ija) experiments were performed in 
a low-turbulence closed wind tunnel of 80 x 80 cm^ cross 
section. The wind velocity U could be varied continu- 
ously up to 65 m s~^ and was measured with a Pitot tube. 
The plates were clamped in a streamlined mast crossing 
the wind tunnel vertically. To ensure that the mast had 
negligible effects on the instability, two masts of differ- 
ent cross section were used (mast 1: thickness x chord== 
4 mm X 20 mm; mast 2: 9 mm x 43 mm). Plates were cut 
in Mylar sheets of mass per unit area m = 140 g m~^ and 
flexural rigidity D = 0.48 x 10"^ N m. To estimate D, the 
deflection due to gravity was measured for horizontally 
clamped Mylar strips of various lengths. 

In all experiments presented in this paper, the same 
protocol has been followed. The plate is clamped in the 
mast and the flow velocity is slowly increased starting 
from zero. At small velocities the plate appears stable, 
i.e. stationary and aligned with the flow. Eventually, 
for a critical flow velocity C/*, the plate flutters sponta- 
neously with a large amplitude and a well-defined fre- 
quency. Then the flow velocity is slowly decreased in 
small decrements until the plate returns to its stable state 
again at the critical velocity C/J as illustrated in Fig. [51 

Visualizations were carried out through the top wall of 
the wind tunnel with a high speed video camera aligned 
with the y-axis (see Fig.[l}i). In the present experiments, 
the camera was operating at 300 Hz with a512x512 pixel 
resolution and the exposure time was varied between 300 
and 1500 fis. The camera visualizations are used to ex- 
tract the mode amplitude A and its angular frequency 
Lu*. This is done by detecting for each snapshot the plate 
deflection at 3/4 of the total plate length. 

The mode amplitude A is plotted in Fig. [2] as a func- 
tion of the reduced velocity and exhibits a hysteretic be- 
havior. This has been observed in all experiments with 
sufhciently large aspect ratio {H* > 1). This hystere- 
sis together with the large mode amplitude observed at 
threshold evokes a subcritical instability bifurcation as 
recently suggested by numerical results llj. However, in 
other nonlinear numerical models 1^, |13| , the bifurcation 




FIG. 2: Mode amplitude A (in arbitrary units) as a function 
of the reduced velocity U* . Amplitude is plotted as filled 
triangles when velocity is increased and open triangles when 
it is decreased. The experimental parameters are H* = 1 and 
M* = 0.6. The instability exhibits a strong hysteresis with 
U* = 10.3 and U^ = 8. However, the amplitudes are well 
fitted by the square-root law A = 58(L^* — 7.4)^" (dashed 
line) suggesting a supercritical bifurcation. 



is supercritical. Figure [2] shows that the amplitudes near 
threshold can be reasonably well fitted with a square-root 
law giving evidence that the bifurcation could indeed be 
supercritical and that the observed hysteresis could just 
be an artifact as discussed below. 

The critical velocities U* and f/J are plotted in Fig. [3] 
as a function of the aspect ratio H* and the mass ra- 
tio M*. First the plate span is varied keeping its length 
constant. This allows to vary the aspect ratio without 
changing the other parameters. As expected from slen- 
der body theory, the critical velocity tends to decrease 
as aspect ratio is increased. From Fig. |3Ua) two regimes 
can be distinguished. For small aspect ratio {H* < 1), 
the measured critical velocities are well predicted by the 
linear stability analysis and the hysteresis is very small. 
For larger aspect ratio {H* > 1) the hysteresis greatly in- 
creases as well as experimental scatter and the linear sta- 
bility analysis tends to underestimate the critical velocity 
U*. Moreover, this threshold seems to slightly increase 
for 1 < H* < 3 contrarily to the decreasing threshold U^ 
and contrarily to the predictions. 

In Fig. ^h) , the aspect ratio has been kept constant 
{H* = 1) by varying by the same factor the plate length 
and span. The measured critical velocity U* is compared 
to the theoretical predictions of the present analysis [10| 
and the 2D model of Kornecki et al. [5[ . In the present 
analysis, the finiteness of the span is taken into account 
leading to a smaller destabilizing term {[p]) in Eq. ^ 
compared to the 2D model and thus a better prediction 
of the instability threshold. Both models allow to pre- 
dict different instability modes as M* is increased as il- 
lustrated by the different lobes in Fig. [3]Jb) . If modes are 
numbered by order of ascending frequencies, the mode 
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FIG. 3: Critical velocities U^ and U^ as a function of the 
aspect ratio H* for M* = 0.6 (a) and as a function of the 
mass ratio M* for H* = 1 (b). Filled squares and triangles 
correspond to measured f/* for the masts 1 and 2 respectively 
and the same open symbols refer to U^. Crosses correspond 
to the experiments of Ref. n fo r 0.6 < H* < 1.5. The solid 
line is the 3D theory of Ref. [lO for same the parameters, the 
dashed line is the slender body theory (similarly to Ref. 0) 
and the dash-dotted line corresponds to the 2D theory with 
circulation (same as Ref. |3 with better computer accuracy). 



two is observed for the smallest mass ratios (Af* < 1.5). 
This single-neck mode is pictured in Fig. IUa,b). For 
larger mass ratios (1.5 < M* < 5), the mode three, a 
double-neck mode, is the first unstable one as illustrated 
in Fig. [Hc,d). Note that the mode one (the mode with 
the smallest frequency) is never unstable for a clamped- 
free plate as it is discussed in refs. [8|,|l0[. For larger mass 
ratios (M* > 5), one would observe higher order modes 
(but in such experiments, plates are so thin that gravity 
makes them twist). In other words, as M* increases, the 
typical mode wavelength A is a smaller fraction of the 
plate length L and therefore if the aspect ratio is kept 
constant, H/X increases. This span to wavelength ra- 
tio has to be large for a 2D theory to be valid and This 
explain why the difference between the 2D and 3D theo- 
retical approaches reduces as M* increases in Fig. ^h). 
In Fig. [ll|a,c), superimposed views of the plate 




FIG. 4: Superimposed views of the flutter modes for two 
plates of same aspect ratio H* = 1 but with different sizes 
[M* = 0.74 for (a) and (b), M* = 1.94 for (c) and (d)]. These 
amplitudes are obtained at threshold and the measured fre- 
quencies are ui* = 1.82 for [/* = 8.1 (a) and lo* = 2.5 for 
f/* = 10.9 (c). The experimental snapshots are compared 
with the modes predicted by the present linear stability anal- 
ysis (with arbitrary amplitude): a;* = 2.4 for V^ = 7.2 (b) 
and w* = 3.7 for [/* = 10.4 (d). 



uniform along the span and exhibit a maximum at mid- 
span and zeros at the plate edges. Second, the gravity 
field induces a non trivial stress tensor in the plate that 
induces compression at the top corner and tension at the 
bottom corner. Since tension has a stabilizing efl^ect, this 
may explain why the top corner is the most unstable po- 
sition. Finally, one should note that bending along the 
span acts as a stiffening effect (because Gauss curvature 
is energetically costly) and therefore it has a stabilizing 
effect on the flutter instability. 

As explained above, the leading-edge pressure singu- 
larity is not properly taken into account in the present 
analysis (mostly because it greatly simplifies the prob- 
lem). This explains why the solid curve in Fig.[3l^a) does 
not tend for large H* to the 2D limit represented by the 
dash-dotted line. If a proper full three-dimensional sta- 
bility analysis was carried out, the critical curve U*{H*) 
would be slightly above the solid line of Fig.[3lja) in order 
to converge to the 2D limit. In this case, the theoretical 
threshold would be even closer to the open symbols rep- 
resenting the descending threshold U^. This is another 
argument for a supercritical bifurcation. Indeed, the in- 
stability could be delayed and the measured U* higher 
than U^ because of small imperfections of the flow (due 
to the wind tunnel or the wake of the mast) driving a 
small two-dimensional plate vibration or deflection and 
thus stiffening the system as exposed above. 

The authors acknowledge support from the French 
ANR (No. ANR-06-JCJC-0087). 



recorded with the camera during one flutter period are 
shown. These visualizations were carried out at the in- 
stability threshold U* for plates of same aspect ratio 
H* = 1 but for two different mass ratios. As predicted by 
the analysis, different modes are observed as M* varies. 
From these data, the mode shapes and their frequencies 
can be compared to the theoretical predictions for the 
same parameters. However nonlinear effects are prob- 
ably important given the large amplitude observed at 
threshold and the agreement with the results of the linear 
stability analysis are only qualitative. 

In the experiments with large aspect ratio H* > 1, 
the flutter mode is no longer purely one-dimensional as 
the top corner flutters with a larger amplitude than the 
rest of the plate. Moreover, in some experiments with 
a very large aspect ratio H* = 4 and for flow velocity 
just below threshold, the same top corner flutters while 
the rest of the plate is still motionless. This indicates 
that the Euler-BernouUi beam equation for the plate is 
no longer valid and one should use the nonlinear Foppl- 
von Karman equations as soon as the deflection along 
the span is of the order of the plate thickness. The rea- 
sons for the two-dimensionality of the plate deflection 
are twofold. First, the pressure field over the plate is non 
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